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Non-statistical event-by-event fluctuations are considered an important signal for the critical end- 
point of the QCD phase diagram. Event-by-event fluctuations of different observables are thus 
investigated in detail in current experiments. In this contribution, a study of the centrality depen- 
dence of event-by-event fluctuations of particle yield ratios by the NA49 experiment in Pb+Pb 
collisions at 158A GeV beam energy is presented for the first time. An increase of the absolute 
values of dynamical fluctuations towards lower centralities is observed. The influence of reso- 
nance decays will be discussed. 

Event-by-event fluctuations are an important observable to be studied at the future CBM experi- 
ment at FAIR. CBM will investigate the intermediate region of the QCD phase diagram in great 
detail searching for the first order phase transition line and the expected critical endpoint. It is 
therefore important to investigate the sensitivity of the CBM detector to particle ratio fluctuations 
in Au+Au collisions at 10-45 AGeV beam energy. Detailed simulation studies will be presented. 
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1. Introduction 

Recent lattice QCD calculations expect a first order phase transition from hadronic to partonic 
degrees of freedom at finite temperature and baryon chemical potential ||]. This first order tran- 
sition line ends with a critical endpoint |^]. The search for either the first order phase transition or 
the critical endpoint is a challenging task in modern high energy heavy-ion physics. 
Due to density fluctuations at the critical point or in the coexistance region at a first order phase tran- 
sition, fluctuations in particle yields and kinematic properties may occur. Indeed, in lattice QCD 
calculations the quark number susceptibility develops a peak at the critical temperature when in- 
creasing the quark chemical potential and thus approaching the expected critical point (figure [T]) [Q]. 
Fluctuations as a measure of susceptibilities [||] should than be enhanced. 




Figure 1: The quark number susceptibility as a function of temperature for different quark chemical poten- 
tials [[§. 

In this work we will concentrate on event-by-event fluctuations of particle yield ratios, such 
as kaon to pion (KJn), proton to pion (p/7T) and kaon to proton (K/p) ratios. Recently, the NA49 
and STAR collaborations have published their results on the energy dependence of the dynamical 
fluctuations of KJti and p/Ti ratios [|| 0J. The measured dependence is illustrated in figure ||. 

The dynamical fluctuations of the KJti ratio increase towards lower collision energies and are 
not reproduced by the UrQMD model, while the p/n ratio fluctuations from the data analysis are in 
good agreement with model calculations. 

In [|]] it was discussed, that this increase of the Kl% ratio fluctuations towards lower beam energies 
could be associated with a decrease of particle multiplicities. Thus one would like to study the cen- 
trality dependence of the fluctuation signal in order to understand the behaviour of this observable 
with the changing average multiplicity of the accepted tracks. In addition, it is neccessary to re- 
measure the low energy region with better experimental resolution and larger acceptance to provide 
high quality data on this observable, which is one of the goals of the future CBM experiment. 
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Figure 2: Energy dependence of the dynamical fluctuations of the ¥Un and p/n ratios in central A + A 
collisions published by the NA49 and STAR experiments 



2. Centrality Dependence of Particle Ratio Fluctuations in NA49 

NA49 is a fixed target experiment at the CERN SPS. Measurements of p + p, p + A and A + A 
collisions were recordered at different energies. NA49 is a large acceptance hadron spectrometer 
which includes two dipole magnets for momentum determination, two vertex TPCs for the mea- 
surement of vertices and particle identification (PID) using dE / dx, two main TPCs for PID with 
dE I dx, two TOF walls for PID at midrapidity using time-of-flight and a Zero Degree Calorimeter 
for centrality determination. A detailed description of the NA49 detector can be found in [JSj] . 
In the current analysis, the specific energy loss dE /dx, measured by the main TPCs of NA49 was 
used for particle identification. The NA49 spectrometer works in the relativistic rise region of the 
Bethe Bloch parametrisation of dE /dx This implies some restrictions on the PID performance, 
namely a lower momentum cut of 3 GeV/c and no possibility of PID on the track4jy-track level. 
Thus, particles have to be identified on a statistical basis by extracting the relative yields in a single 
event using the Maximum Likelihood Method. The data analysis technique is described in detail 



in [lC.y]. 



In this section we will focus on the measured centrality dependence of the dynamical fluctuations 
of particle yield ratios. The centrality of an A + A collision was determined by measuring the 
total energy of the projectile spectators deposited in the Zero Degree Calorimeter of NA49. The 
dynamical fluctuations of a particle yield ratio are defined as the geometrical difference between 
the relative width (a=RMS/MEAN) of the eventwise ratio distribution for data (Odata) and mixed 
(<W) events: 
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(2.1) 



2.1 Centrality Bin Size 



The main motivation to use particle yield ratios in the fluctuation analysis is that in first ap- 
proximation volume fluctuations are cancelled in the ratio. Nevertheless, the fluctuations of the 
projectile and target participants in Pb + Pb collisions may have an influence on the measured sig- 
nal, thus more detailed studies had been performed on the subject of the dependence of dynamical 
fluctuations on the selected centrality bin size. In the published NA49 results on the energy depen- 
dence of the 3.5% most central Pb + Pb collisions were used in the analysis [||, [T0| ]. 
The distribution of the total energy of the projectile spectators deposited in the Zero Degree Calori- 
meter of NA49 is shown in figure || (left plot). The centrality was integrated in the following bins: 
(0-3)%, (0-3.5)%, (0-5)%, (0-20)% most central Pb + Pb collisions at 158A GeV beam energy. 
The observed dependence of the K/7T ratio fluctuations is shown in figure |5| (right plot). 
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Figure 3: Left picture: distribution of the energy deposited in the Zero Degree Calorimeter (ZDC) of NA49 
in Pb + Pb collisions at 158A GeV beam energy. Vertical lines show cuts, which correspond to 3%, 3.5%, 
5%, 10%, 15%, 17.5% and 20% centtality bin size. Right picture: dynamical fluctuations of the KJn ratio 
as a function of the centrality bin size in Pb + Pb collisions at 158A GeV beam energy. 



The results are in good agreement with the top SPS beam energy point from [^]. As has been 
discussed before, the fluctuations of the number of projectile and target participants have an effect 
on the measured signal, and a slight increase of about 2% of fluctuations with increasing integrated 
centrality of Pb + Pb collisions is observed as a general trend. However, the difference between the 
values of the dynamical fluctuations of the kaon to pion yield ratio in case of using 3.5% and 5% 
most central Pb + Pb collisions is small. So, in order to gain statistics, which is important for the 
errors, the centrality dependence of the dynamical fluctuations can be studied using a centrality bin 
size of 5%. 
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2.2 Simulations with UrQMD 

Particle identification on the statistical basis can introduce a bias to the event-by-event fluctua- 
tions. This bias might be larger in semi-peripheral Pb + Pb collisions, since in this case the average 
track multiplicity decreases and the identification procedure becomes more difficult. It becomes 
particularly challenging if single particle multiplicities in the acceptance approach zero. As no 
negative numbers can be allowed for stable fit results, the eventwise KJn distribution develops a 
spike at zero [Oj ]. 

In order to find out the working limits for the analysis of the centrality dependence of Pb + Pb 



collisions, simulations with the UrQMD model [12] have been performed. The centrality of a col 



lision was determined using the provided value of the impact parameter. Particles, generated by 
the model, were processed through the acceptance filter of NA49, which incorporates the specific 
track cuts used in the current data analysis. For each accepted track, a dE/dx value was simulated 
using a parametrization of the measured dE / dx distribution in NA49. The accepted tracks with 
their simulated dE/dx response were processed through the same analysis routines as the tracks 
from data, namely the event-by-event fit of the eventwise dE/dx distribution (E-b-e fit). As an un- 
biased reference, particle identification based on the particle type provided by the UrQMD model 
was used, so called Monte Carlo identification (MC PID). The extracted centrality dependence of 
the dynamical fluctuations of the kaon to pion yield ratio for these two cases is shown in figure ^. 
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Figure 4: Dynamical fluctuations of the K/7T ratio as a function of the number of wounded nucleons in Pb + 
Pb collisions at 158A GeV beam energy, simulated with UrQMD, for two cases: Monte Carlo based counting 
(triangles) and event-by-event fit based on simulated dE / dx response (squares). 



For both identification procedures the dynamical fluctuations increase with decreasing num- 
ber of wounded nucleons. But the results for the event-by-event fit start to deviate from the MC 
reference for less than 200 wounded nucleons, which approximately corresponds to a centrality of 
35%. The contribution of the spike at zero in this case exceeds 2% of all measured events. Thus 
NA49 data will be analyzed in the centrality range from 0% to 35% most central Pb + Pb collisions, 
where the identification procedure does not introduce a bias to the results. 
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2.3 Centrality Dependence 

The measured centrality dependences of the dynamical fluctuations of the KJn and p/n ratios 
in Pb + Pb collisions at 158A GeV beam energy are presented in figure ||. 
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Figure 5: Dynamical fluctuations of the KJn (left picture) and the p/n (right picture) ratios as a function of 
the number of wounded nucleons in Pb + Pb collisions at 158A GeV beam energy. The shaded band shows 
systematic errors determined from a variation of ttack cuts. 

Figure ^ shows the dynamical fluctuations of the K/p ratio as a function of the number of 
wounded nucleons in Pb + Pb collisions at 158A GeV beam energy. The fluctuations are calculated 
as mean value from analysis with different track cuts. Their difference is used to estimate the 
systematic error which is given by the shaded band. 
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Figure 6: Dynamical fluctuations of the K/p ratio as a function of number of wounded nucleons in Pb + Pb 
collisions at 158A GeV beam energy. The shaded band shows systematic errors determined from a variation 
of track cuts. 

The general trend is the same for all considered particle yield ratio fluctuations: the absolute 
value of the dynamical fluctuations increases towards lower centrality of Pb + Pb collisions. 
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2.4 Scaling of Dynamical Fluctuations 



From the definition of the dynamical fluctuations (equation 2.1 ) one can derive an exact ana 



lytical expression, which in case of KJk consists of kaon, pion and correlation terms: 



var{N K ) var{N n ) cov{N K ,N n ) 
° dyn = V " 2 <N K XN X >> M 

In case of p/n, Nk should be substituted with N p . As has been argued in ||6[ [13|], dynamical 
fluctuations of the pin ratio originate dominantly from A resonance decays, which produces corre- 
lated protons and pions. The analytical expression for the dynamical fluctuations in this case can 
be approximated by the correlation term only: 



y/n^ cov{N p ,N n ) = {<N p ><N n >r 

d y n ~ ]j<N p ><N x > V <N p ><N n > ' K ' J 

where we have parametrised the covariance by the product of average multiplicities of protons 
and pions to the power of a. With a strong feeddown from resonances this parameter a is expected 
to become 0.5. In case of the KJn ratio we have Nk <C N n and the analytical expression for Od yn 
can be apporoximated by: 



K/ n var(N K ) 
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Figure 7: Left picture: dynamical fluctuations of the KJn ratio as a function of < Ng > for the energy 
(circles) and the centrality (squares) dependences. Right picture: dynamical fluctuations of the p/n ratio as 
a function of < N p >< N K > for the energy (circles) and the centrality (squares) dependences. Solid curves 
show fits according to equation 2.3. 

Figure [7] shows the dynamical fluctuations of the KJn ratio as a function of < Nk > (left pic- 
ture) and the dynamical fluctuations of the p/n ratio as a function of <N P ><N 7[ > (right picture). 
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Data points are taken from the measured energy and centrality dependences for a more loose set of 
track cuts; the mean number of particles is calculated in the acceptance used to determine the fluc- 
tuations. After fitting the dependences of the pin ratio fluctuations with equation [2.3| , the following 
a parameters were obtained: a = 0.66 ±0.12 for the energy dependence and a = 0.51 ±0.03 
for the centrality dependence. This observation strongly supports the assumption that A resonance 
decays are the dominant source of dynamical fluctuations of the p/n ratio and give a natural expla- 
nation for both, the energy and centrality dependence. As expected also for KJn ratio fluctuations 
some dependence on < Nk > is seen, however the energy and centrality dependence obviously 
scale differently with < Nk >■ Whether this might be connected to different relative strangeness 



production in smaller systems [14] requires further studies. 



3. Particle Ratio Fluctuations in CBM 



CBM is the future heavy-ion fixed target experiment at FAIR []15[]. Its goal is to measure 
hadrons and leptons in A + A collisions at beam energies from 10 - 45 AGeV. The proposed CBM 
detector (electron option) is shown in figure |8[ 




Figure 8: Sketch of the CBM detector. 

In this section we will focus on feasibility studies of the measurement of event-by-event fluc- 
tuations of the KJn ratio with CBM. In CBM hadrons will be identified using a TOF wall. With 
a time resolution of better than 80 ps it will be possible to identify particles on the track-by-track 
level, which together with a large acceptance (larger as compared to the NA49 experiment) should 
avoid the development of a spike at zero in the eventwise KJk ratio distributions. The distribu- 



tion of the squared mass versus momentum of reconstructed primary hadrons from UrQMD [ |12| ] 
simulations is shown in figure ^. 

A clean identification of kaons with a pion contamination on the level of 15% only (integrated) 
is possible up to momenta of approximately 3 GeV/c. 

A certain level of the purity of kaon identification is achieved by using an upper momentum 
cut. As this of course restricts the phase space acceptance, the effect of this cut on the dynamical 
fluctuations of the KJn ratio was investigated. As a reference, the results with Monte Carlo (MC) 
identification but the same upper momentum cut were calculated. The observed dependence of 
the dynamical fluctuations on the purity of kaon identification i.e. the applied momentum cut is 
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Figure 9: Distribution of the squared mass versus momentum of reconstructed primary hadrons from central 
Au + Au collisions at 25A GeV beam energy (UrQMD simulation). Right picture shows a projection for the 
momentum of 3 GeV/c. The tails originate from secondary particles and ghost tracks (mismatches). 

presented in figure [T0| for MC identification and fully reconstructed and identified tracks in the 
CBM simulation. 
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Figure 10: Dynamical fluctuations of the K/n ratio as a function of purity of the kaon identification for 
the standard identification procedure (circles), Monte Carlo based counting with the upper momentum cut, 
which corresponds to the chosen purity, (squares) and full CBM acceptance (filled region). 



For example the purity requirement of at least 50% per momentum bin corresponds to a mo- 
mentum cut of 6 GeV/c, 99% to 3.25 GeV/c. Values of the dynamical fluctuations increase if the 
momentum cut is more restrictive. Good agreement between results with real identification and 
PID on the MC level with corresponding momentum cut is observed. The results are closer to the 
values for the full CBM acceptance in case of lower purity, which corresponds to a more relaxed 
total momentum cut. 
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4. Summary and Outlook 

The centrality dependences of KJn, pin and KJp ratio fluctuations show increasing absolute 
values towards lower centrality. A multiplicity dependence is indeed expected and the dependence 
of the p/n ratio fluctuations on energy and centrality can both together be well described assuming a 
dominant contribution from A resonance decays. In case of the KJn ratio fluctuations however, the 
values of the dynamical fluctuations of the KJn ratio for semi-peripheral Pb + Pb collisions at 158A 
GeV are higher than those for central Pb + Pb collisions at lower SPS energies at the same values 
of < Nk >■ In other words: the energy and centrality dependences of the KJn ratio fluctuations 
scale differently with <Nk>- 

As dynamical fluctuations in particle emmision are expected to be an important observable for 
the critical endpoint and current data are not conclusive yet, the future CBM experiment at FAIR 
prepares to also provide precise evaluation/observation of these fluctuations. Feasibility studies for 
the measurement of event-by-event fluctuations of the KJn ratio in central Au + Au collisions were 
presented and demonststrate that some of the experimental difficulties in the NA49 analysis can be 
overcome. The study shows that such a measurement is feasible and that no strong bias from the 
identification procedure with the TOF wall of CBM is expected. 
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